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Abstract 

Galactic halo gas traces inflowing star formation fuel and feedback from a galaxy's disk 
and is therefore crucial to our understanding of galaxy evolution. In this review, we 
summarize the multi-wavelength observational properties and origin models of Galactic 
and low redshift spiral galaxy halo gas. Galactic halos contain multiphase gas flows that 
are dominated in mass by the ionized component and extend to large radii. The densest, 
coldest halo gas observed in neutral hydrogen (HI) is generally closest to the disk (< 20 
kpc), and absorption line results indicate warm and warm- hot diffuse halo gas is present 
throughout a galaxy's halo. The hot halo gas detected is not a significant fraction of a 
galaxy's baryons. The disk-halo interface is where the multiphase flows are integrated 
into the star forming disk, and there is evidence for both feedback and fueling at this 
interface from the temperature and kinematic gradient of the gas and HI structures. 
The origin and fate of halo gas is considered in the context of cosmological and idealized 
local simulations. Accretion along cosmic filaments occurs in both a hot (> 10 5 ' 5 K) 
and cold mode in simulations, with the compressed material close to the disk the coldest 
and densest, in agreement with observations. There is evidence in halo gas observations 
for radiative and mechanical feedback mechanisms, including escaping photons from the 
disk, supernova-driven winds, and a galactic fountain. Satellite accretion also leaves 
behind abundant halo gas. This satellite gas interacts with the existing halo medium, 
and much of this gas will become part of the diffuse halo before it can reach the disk. 
The accretion rate from cold and warm halo gas is generally below a galaxy disk's star 
formation rate, but gas at the disk-halo interface and stellar feedback may be important 
additional fuel sources. 
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1 Introduction 



Halo gas connects the baryon-rich intergalactic medium (IGM), to the star- 
forming disks of galaxies. It represents a galaxy's future star formation fuel 
and the result of galactic feedback processes. A galaxy's halo gas is therefore a 
combination of new material from the IGM and satellite galaxies and recycled 
material from the disk. The metallicities of the long-lived stars in the Galactic 
disk indicate the vast majority of the incoming fuel should be low metallicity, and 



therefore feedback material cannot be the dominant fuel at all redshifts ( 


Chiap- 


pini et al., 2001 


Larson, 1972 


Fenner & Gibson, 2003; 


Sommer-Larsen et al., 



2003). In addition, star formation rates (SFRs) indicate that ongoing accretion 



is required for galax ies at a range of redshifts (Erb, 2008 Putman et al., 2009c 
Hopkins et al., 2008); the Milky Way is no exception with only ~ 5 x 10 9 Mq of 



fuel in the disk and a current SFR of 1-3 M Q /yr (Chomiuk & Povich, 2011 



Smith et al., 1978; Robitaille & Whitney, 2010). All new galactic fuel needs to 



pass through a galaxy's halo, where the gravitational potential will draw the gas 
to the disk and help maintain spiral structure by keeping the Toomre Q parameter 
low ( |Toomre, 1990] |1964[ ). 

This review outlines the observational properties and potential origins of the 
multiphase halo medium in both the Milky Way and other spiral galaxy halos. 
Though a multiphase halo was first introduced theoretically by Spitzer (1956), 
the observational relationship between the multiple phases has only become more 
apparent in the past decade. Low redshift studies of halo gas are key, as this is 
where physical details can be resolved and the low density gas can be mapped in 
emission. We define halo gas to lie within the approximate virial radius (~ 250 
kpc for the Milky Way) and beyond the star forming disk. The latter boundary 
is often difficult to define and most likely somewhat artificial given the ongoing 
transition of material between the lower halo and the disk. This disk-halo inter- 
face is discussed at several points in the review since it represents a key transition 
point. The total accretion rate at this interface is also an important quantity to 
determine, given questions about a galaxy's future star formation and evolution. 

There are several possible origins for halo gas and a mixture of origins is likely 
for a spiral galaxy like the Milky Way at low redshift. In cosmological simulations 
of galaxy formation, a large amount of the halo gas originates from inflowing IGM 
along cosmic filaments. These simulations also find galactic feedback provides 
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enriched material to galaxy halos. Finally, satellites contribute gas to galaxy 
halos as they accrete onto the galaxy and are stripped of their gas. As we are 
now able to determine the physical properties of much of the multi-phase halo 
medium (with distance determinations and detailed observations), origin models 
can be more tightly constrained. In addition, simulations are increasingly able 
to include the relevant gas physics and guide future observations. 

Through observations and simulations we have made substantial progress in 
understanding gaseous halos and the flow of baryons in the universe. We begin 
this review with an overview of observations of the gaseous Galactic halo (Q. 
Due to its proximity, this is the halo for which we have the most data, and we 
discuss its gas in sub-sections that are largely defined by temperature. The gas 



at the disk-halo interface of the Milky Way is discussed separately in £ 2.5 The 
Milky Way has a total mass of ~ 2 x 10 12 M Q (|Reid et al., 2009| |Kalberla et al., 



2007), and ideally we would discuss only galaxies of this mass for comparison. In 
reality, there is no perfect nearby Milky Way equivalent, or limited observations 
of the best examples, and so ^3] summarizes the halo gas properties of a range 
of low redshift spiral galaxies. We consider the origin of the observed halo gas 
in the context of simulations and theoretical models in £j4] and comment on the 
survival of halo clouds in §[5j Finally, the sources available and possible methods 
of feeding galaxy disks are presented in £j6] and the outlook for the future is noted 
in fj7) 

2 The Galactic Halo 

This section focuses on the significant observational progress that has been made 



in studies of Milky Way halo gas in the past decade (see Wakker & van Woerden 



(1997) for earlier work) and is divided into sub-sections based primarily on the 
temperature of the gas. The neutral hydrogen halo gas discussed in §2.1| is at tem- 
peratures < 10, 000 K and detected through the 21-cm line in emission. This is 
traditionally the primary tracer of Galactic halo gas, and we therefore separately 
discuss its physical properties (§ |2.1.1 ). The warm and warm- hot gas is defined 



here as largely ionized with temperatures > 10 4 K, but < 10 6 K (^2.2). This gas 
is detected in optical line emission (primarily Ha) and absorption lines in the 
optical and ultraviolet . Finally, hot gas is at temperatures greater than 10 6 K 
and detectable through x-ray observations close to the disk and indirect methods 
further out (32.31). The dust, molecules and metals in halo gas are discussed in 



52.4 



Distances are often an issue when defining g halo medium. There has 



been significant progress in this area (see £2.1.1), but in any case, it should be 
noted that halo gas has largely been defined based on its velocity not falling 
within the expected Galactic disk velocity in a given direction. This definition 
may in some (relatively rare) cases pick up Galactic gas that has been recently 
imparted with a large amount of kinetic energy, and it definitely misses halo gas 
that overlaps with the velocity of the disk. In general, the transition between 
halo and disk gas is blurry, and a continuum of gas between the disk and halo 
seems likely. Halo gas is defined here to be at heights above the Galactic Plane 
(z-heights) of a few kpc and not part of the Galactic warp, and gas below this 
that shows some connection to the disk is discussed in the disk-halo interface 
section (pl>|). 
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2.1 Neutral Hydrogen Halo Gas 



The first halo gas detected was in neutral hydrogen emission at high velocity 
relative to Galactic disk velocities (Muller et al., 1963). These discrete gas clouds 
were subsequently referred to as high- velocity clouds (HVCs). HVCs are the 
coldest and densest component of halo gas and are found throughout the sky 
(Figures [l] & [2]) . They are material stripped from satellites and/or cooling com- 
ponents of the accreting IGM and feedback material (see Q . Significant progress 
in mapping the distribution of HVCs has come through multibeam instruments 
that have been used to observe large areas of sky at relatively high resolution 
(4' to 16'; [Putman et al., 2002[ |Saul et al., 2012[ |McClure-Griffiths et al., 2009| ). 
In addition, the entire sky has now been mapped at a uniform reso lution of 
35' and 1 km s _1 with the Leiden/ Argentine/Bonn (LAB) HI survey (Kalberla 
et al., 2005, ; see Figures [T] & [2]) . The full-sky successor to the LAB survey will 
be a combination of multibeam surveys with a resolution of 9 — 16' and 1 km 
s _1 (iMcClure-Griffiths et al. 



2009 



Winkel et al., 2010[ ). 
HVCs range in size from 1000 deg 2 to the resolution limit of the observations, 
and large clouds are often found to be composed of a number of smaller clouds 



when observed with better resolution and equivalent sensitivity (Putman et al., 



2002; Stanimirovic et al., 2008 Hsu et al., 2011). There are many catalogs of 



clouds that have been created with the various HI survey data that use different 
catalog methods. The number of individual clouds in these catalogs ranges from 



600 over the entire sky in the lowest resolution surveys (35' beam; Wakker 
fc van Woerden, 1991 ), to 2000 in just the southern sky (15.5' beam; |Putman" 
et al., 2002[), to ~ 700 clouds smaller than 20' in 7600 deg 2 of the sky (4' beam; 



Saul et al., 2012). It remains difficult to separate the Galactic disk gas from halo 
gas. One method of cataloguing distinct halo features is to run algorithms on the 
entire HI data cube to determine structures that are distinct from each other in 



position- velocity space (Mailer et al. in prep.; Sharma & Johnston, 2009). 



Clouds of all sizes can be grouped into complexes of most likely related clouds 
based on their spatial and kinematic proximity as labeled in Figures [T] and [2] In 
the Southern hemisphere the HVCs that dominate in sky coverage and mass are 
associated with the Magellanic System. In particular, the complexes labeled MS 
and LA are the trailing Magellanic Stream and Leading Arm (recent references 



include: Putman, 2004[ |Putman et al., 2003bl 


1998| IBriins et al., 2005; 


Wcst- 


meier & Koribalski, 2008; iStanimirovic et al., 2002 


|2008| [Nidever et al., 2010 


)• 



These complexes were created from the interaction of the Large and Small Magel- 
lanic Clouds (LMC and SMC) with each other and with the Milky Way. The 
Magellanic Stream is a long, continuous structure with a well-defined velocity and 
column density gradient, while the Leading Arm is a collection of clouds through- 
out the region leading the Magellanic Clouds. In the Northern hemisphere there 
are a number of complexes of similar spatial size, but Complex C is the largest. 
The best way to view the HI component of the complexes at their approximate 
distances (see 2.1.1) is through the three-dimensional Figure [3j In both hemi- 



spheres there are a number of small clouds called compact HVCs (CHVCs) that 
have been claimed to be at large distance due to their size (< 2° in diameter) 
and relative isolation ( Braun &: Burton, 1999} |de Heij et al., 200*2 ) . CHVCs can 



largely be associated with known HVC complexes in position-velocity space as 
shown by the plus symbols in Figures [T] and [2} and are therefore unlikely to be at 
greater distance. Even the ultra-compact HVCs (< 20') have been shown to be 
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related to HVC complexes (ISaul et al., 2012h. None of the CHVCs or HVC com- 



plexes have a clear association with stellar features in the halo (Hopp et al., 2007 



Siegel et al., 2005), excluding the association of MS and LA to the Magellanic 



Clouds. 

The majority of the kinematic range of the HVCs can be seen in Figures [T] 
and [2] The range has been extended further with new surveys and is currently 
approximately —500 to +450 km s -1 in the Local Standard of Rest velocity 
frame (Vlsr) or between —300 and +300 km s _1 in the Galactic Standard of 
Rest frame (Vgsr)- HVCs can be defined by the cloud's deviation velocity (V(j ev ), 
or the amount the velocity of the cloud deviates from a simple model of Galactic 
rotation. Figur es[l] and [2] were created with |Vd cv I > 75 kms -1 and thus there are 
some small differe nces from maps created with the classical | Vlsr I > 70 — 90 km 
s _1 definition (e.g. |Westmeier, 2007 Wakker & van Woerden, 1991). The clouds 
at both the positive and negative velocity extremes are largely associated with 
the Leading Arm and the tip of the Magellanic Stream, respectively. Since HVCs 
have not been catalogued in the velocity range obscured by Galactic emission, it 
is difficult to assess the dispersion of the population. Using the GSR frame, an 
estimate for the velocity dispersion of the HVCs is ~100-120 km s _1 . 

The HVCs have typical linewidths of 20-30 km s -1 , corresponding to velocity 
dispersions of a = 8.5 — 13 km s _1 ( Putman et al., 2002 de Heij et al., 2002 



Kalberla & Haud, 2006 


. This linewidth 


for HVCs at ~ 9000 K 


Hsu et al., 2011 



Many HVCs also show a n arrow com- 
ponent in their velocity profiles corresponding to cold gas at < 500 K (Kalberla & 



Haud, 2006). The two-phase structure of HVCs is expected for clouds in pressure 



equilibrium with a diffuse hot Milky Way halo medium (Wolfire et al., 1995:2.3 



It is possible that HVCs that do not show this two-phase structure are currently 
in a more turbulent environment, or are embedded in a lower pressure part of 
the halo. 

HVCs have typical peak column densities of ~ 10 19 cm -2 , and a small number 
of clouds have peak column densities above 10 20 cm -2 with the 4-15.5' resolu- 
tion HI surveys (Putman et al., 2002 Hsu et al., 2011). Additional dense cores 

10 2U — 



cm have been detected within some HVCs 



Wakker et al., 2002; Ben Bekhti et al., 2006). 



with HI column densities > 
with synthesis observations (e.g., 
The clouds are found to extend down in column density to the limits of the HI 
observations (> 10 17 cm~ 2 ; |Lockman et al., 2002HBraun & Thilker, 2004[). This 



suggests HI surveys may be detecting only a fraction of the area of many com- 
plexes, although they are likely to be detecting the majority of the HI mass (e.g., 



Peek, 2009). There is also evidence that the clouds detected in HI emission are 
linked to lower column density absorption line systems (Ben Bekhti et al. 
Murphy et al., 19951 iNidever et al., 20101). At ~ 7 x 10 17 



2012 



7 x 



cm in HI emission 



the covering fraction is ~ 37% of the sky (Murphy et al., 1995 Lockman et al., 



2002). 

When the detailed HI structure of HVCs is examined, the clouds show evid- 



ence for being ionized by an external radiation field (Maloney Sz Putman, 2003 



Sternberg et al. , 2002 ) . The HI structure also often indicates the clouds are mov- 



ing through a diffuse halo medium and being disrupted (Peek et al., 2007 Brims 
et al., 2000; Stanimirovic et al., 2006, see Figure[4]and ^5). Head-tail clouds, that 



form as the leading side of the HVC is compressed and gas is stripped behind 



it, are particularly abundant in the vicinity of the Leading Arm (Putman et al., 
201 lbp. 
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2.1.1 Distances and Physical Properties HVC distance measurements 
allow for estimates of their physical properties such as mass, size, volume density, 
and pressure 



(e.g. 



Wakker et al., 2001; Hsu et al., 2011). The best method to 



measure the distance to a HVC is to use a halo star of known distance toward the 
cloud and search for absorption lines at the HVC's velocity in the star's spectrum. 
A detection of absorption from the HVC shows that the cloud is in front of the 
star, while not detecting absorption may indicate the HVC is behind the star. 
Usually the optical Ca II H and K and/or Na II doublet absorption lines from 
the HVC are used for this distance measurement. Concluding the HVC is behind 
the star with a non-detection requires very deep observations. This is because 
the line of sight to the star may have a lower HI column density than measured 
by the radio observations and there is uncertainty in the gas-phase abundance 
of the metals (Wakker et al., 2001). The identification of halo stars with the 
Sloan Digital Sky Survey (SDSS) and other surveys has enabled distances to be 
obtained for almost all of the large, labeled HVC complexes in Figure [TJ Most of 
the complexes are at distances between 2-15 kpc and z-heights ^ 10 kpc (Wakker 



et al., 2008 


Thorn et al., 2006, 


2008, 


Wakker et al., 2007 



and, as viewed from the Galactic Center, they are within 30° of the disk plane. 
The clear exceptions to these distances are the Magellanic Stream and Leading 
Arm, which are at ~ 55 kpc near the Magellanic Clouds and may extend to 100- 
150 kpc at the tip of the Stream (e.g., Besla et al., 2010). The complexes can be 
viewed at their distances in the three-dimensional map of Figure |3| 

Distance measurements for HVCs obtained with halo stars are referred to as 
direct distance constraints. Indirect distance constraints are model dependent 
and available for almost all of the HVCs. Ha observations of HVCs can be used 
to estimate their distance if one assumes that the Ha emission is the result of 
the ionizing radiation from the Galaxy reaching the surface of the cloud (the 
fraction escaping normal to the disk is f esc ~ 6% in this model; Bland-Hawthorn] 



fc Maloney, 1999 2002 Putman et al, 2003a Bland-Hawthorn fc Putman, 20011 



Bland-Hawthorn et al., 1998). Most of the distances derived from Ha observations 



of HVCs have been consistent with their direct distance constraints using halo 
stars. Clouds without direct distance constraints (in particular the CHVCs) are 
also usually detected in Ha emission, which is consistent with them being within 
40 kpc of the Galactic disk ( |Tufte et al., 2002||Putman et al., 2003a[ ). The striking 
exception to the escaping radiation model is the Magellanic Stream, which is too 
bright in Ha emission given its distance (see {2.2). 

The most important additional indirect constraint on HVC distances is deep 
HI observations of systems similar to the Milky Way and/or Local Group. If 
one assumes the distribution of HVCs is similar, searches of Local Group analogs 
for HI clouds place even the compact HVCs within 80 kpc of the Galactic disk 
(Pisano et al., 2007) and deep observations of Andromeda suggest the clouds 
are within 50 kpc (Westmeier et al., 2007). The distance constraints have ruled 
out the possibility that the majority of the small HVCs are dark matter halos 
found throughout the Local Group (Blitz et al., 1999 Braun & Burton, 1999). 
Recent Arecibo surveys are finding smaller clouds that do not yet have distance 
constraints (Giovanelli et al., 2010 Begum et al., 2010), but most of these clouds 
are associated with HVC complexes in position- velocity space (Saul et al., 2012). 

The HI masses and sizes of HVC complexes can be determined once distances 
are known. Calculating the HI mass requires a measurement of the total flux 
of the cloud and is proportional to the square of the distance. Determining 
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the physical size of a HVC requires scaling the observed size directly with the 
distance, and this measurement may vary with the beam size and sensitivity of 
the observations. The HI masses of the HVC complexes (each consisting of many 
individual clouds) are in the range of 10 5 — 5 x 10 6 M , and their sizes are a few to 
15 kpc across (excluding MS and LA). There are now enough distance constraints 
that we can determine the total mass in halo clouds. For the 13 complexes with 
direct distance constraints, we calculate the mass by placing each complex at the 
average of its upper and lower distance limits. We place the remaining clouds 
(~ 23% of the total HVC HI flux) at a distance of d=10 kpc, since the majority of 
the complexes are at approximately this distance and indirect distance constraints 
indicate there are not clouds at significantly larger distances. This gives a total 
HI mass in non-Magellanic System HVCs of 2.6 x 10 7 M . This mass should 
b e ro ughly doubled and multiplied by 1.4 to account for the ionized component 
(£2.2) and Helium respectively, thus giving a total mass of 7.4 x 10 7 Mq. The 
expected ongoing accretion rate from the HVC complexes is addressed in §6.1| 
The Magellanic System alone contributes a minimum of 3 x 10 s M (d/55 kpc) 2 
just in HI to the halo from the Magellanic Stream, Leading Arm and Magellanic 
Bridge which joins the two Clouds (Putman et al., 2003b Brims et al., 2005). 
Including the LMC and SMC galaxies would add another 7 x 10 s M in HI to 
the Milky Way halo. All other satellites with HI gas are beyond the virial radius 
of our Galaxy ( Grcevich &; Putman, 2009 ) . 

The derivations of volume densities and pressures require additional assump- 
tions about the HVCs and therefore should be adopted with some caution. For 
volume densities, it is often assumed the cloud's path length is similar to the size 
of the dimensions observed. For spherical clouds this may be a valid assumption 
if the linewidth does not seem unusually broadened. Volume density estimates 
for clouds with direct distance constraints are ~ 0.05 — 0.15 cm -3 , and when 
simulated halo clouds are 'observed' the measured volume densities are within 



a factor of 3-4 of their actual volume densities (Wakker et al., 2001 Hsu et al., 



2011). An estimate of a HVCs pressure is made with the ideal gas law using the 



derived volume density and the kinetic temperature calculated from the linewidth 
(see {2.1). Using the above range of densities and a temperature of 9000 K typ- 
ical for the HVCs' warm neutral component, the pressures for the clouds with 
distance constraints are log(P/k) = 2.7 — 3.1 K cm -3 . The HVCs are at z-heights 
from 3-9 kpc and the pressures derived are consistent with the calculations by 



Wolfire et al. (1995) for clouds in pressure equilibrium with a hot halo medium. 



The estimates for the pressures associated with the Magellanic Stream clouds 
are substantially lower, and this would be expected for clouds in a more distant, 
lower density halo medium ( Stanimirovic et al., 2008; Hsu et al., 2011). 



2.2 Warm Ionized and Warm-Hot Halo Gas 

The warm gas (~ 10 4-5 K) in the Galactic halo is detected with deep Ha emis- 
sion line observations and low and intermediate ion absorption lines (e.g., Si II, 
Si III, C II, C III, O I) using background QSOs or stars. The Ha observations 
are generally targeted toward known HI HVCs, while the location of the absorp- 
tion line observations are dictated by the positions of background probes. Both 
methods have led to numerous detections of warm Galactic halo gas. This gas 
probes the interaction of HI clouds with photons from the disk and with other 
halo gas ( ^4.2[ ), and also comes directly from the IGM and/or satellites (j ]4.1| fc 
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4.3). 



There are numerous detections of HI halo clouds in Ha emission with Fabry- 



Perot instruments and deep long-slit spectroscopy ( Putman et al., 2003a[ Tufte 
et al., 19981 [2002l|Bland-Hawthorn et al., 19981 IWeirier & Williams, 19961 IWeiner 



et al., 2002). All of the high- velocity Ha detections are located on or near HI 



HVCs. Extensions of the Ha emission from the HI have been noted for one 
complex (Hill et al., 2009), but this does not appear to be a generic condition 
( |Haffner, 2005[ |Haffner et al., 2001[|Putman et al., 2003a| ). The detections of Ha 
emission are typically < 0.5 Rayleigh (R), where 1 R is 10 6 /47r photons cm -2 s _1 
sr _1 . There are some non-detections of Ha emission for HVCs, but it is possible 
they will be detected in the future if the cloud is fully mapped to deep levels. One 
distinct non-detection may be the Leading Arm of the Magellanic System, which 
is not detected in Ha for numerous pointings (< 30—50 mR). This is in contrast to 
the Magellanic Stream, which is particularly bright in many locations and requires 
models that invoke collisional ionization mechanisms beyond photoionization by 



radiation from the disk (Bland-Hawthorn et al., 2007; Weiner Sz Williams, 1996 



Putman et al., 2003a; 2.1.1 and £4.2). The velocities of the HI and Ha lines are 



approximately the same for all HVC complexes, and there is no evidence for a 
correlation between HI column density and Ha brightness. Other species such 



as [Nil] and [SII] have also been detected in emission for some HVCs (Bland- 
Hawthorn et al., 19981 iTufte et al., 19981 iPutman et al., 2003a|). 



The mass of the warm ionized component of halo clouds should be added to 
the neutral hydrogen mass to obtain the total mass in halo clouds. The Ha 
detections have been used to calculate that the ionized component of the HVC 
complexes is on the order of the mass of the neutral component in many cases 
( |Hill et al., 2009| |Wakker et al., 2008[ |Shull et al., 20111) . This suggests the warm 
ionized component of the HVCs contributes 3 — 4 x 10 Y Mq to the Galactic halo 
(excluding the MS and LA; see { 2.1.1). It should be noted that the ionized gas 
masses vary by significant factors depending on if one assumes the ionized and 
neutral gas are mixed together, or that the ionized gas forms a skin around the 
neutral component. The numerous Ha detections along the Magellanic Stream 
indicate this feature also has a significant amount of ionized gas associated with 
it (see also Fox et al., 2010). 

Absorption line observations are able to detect much lower column density gas 
at a range of temperatures compared to emission line observations. An important 
conclusion from the absorption line observations is that multi-phase high-velocity 
gas covers a significant fraction of the sky. A covering fraction of ~ 81% was found 



for T~ 10 4 4,5 K gas using high velocity Si III lines ((Nsi m = 2.6 x 10 13 cm 2 ) 



Shull et al., 2009 



was found for T'- 



Collins et al., 2009), and a similar covering fraction (> 60%) 
10 a " b K gas using O VI absorption (No vi > 2.5 x 10 13 cm 2 ; 



Sembach et al., 2003). These results can be used to infer that low column density 



ionized hydrogen covers the majority of the sky (see also Lehner et al. 2012). 
Models show that the Si III absorbers represent photoionized material, while the 
O VI is largely from collisionally ionized gas produced as cool clouds interact 
with a hot halo medium (Sembach et al., 2003). A combination of photo- and 
collisional ionization is commonly needed in models to reproduce the line ratios 
for a ll sightlines (|Shull et al., 2011| |Fox et al. , 20061 |2010l |Haffner & Madsen 



2011). Lehner & Howk (2011) recently constrained the distance to a collection 



of ionized HVCs (IHVCs) using various ultraviolet absorption lines detected in 
the spectra of halo stars. They found the ionized clouds with |Vl5_r| < 170 km 
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s -1 are within ~ 15 kpc and contribute on the order of 10 8 Mq in gas to the 
halo (assuming a metallicity of 0.2 solar; see also Shull et al. 2009). Thus, the 
largely ionized component of halo gas most likely dominates in mass over the HI 
component, although one must be careful that the same velocity cuts are made 
with these comparisons. 

In Figures [T] and [2j high- velocity absorption lines are shown over the HI map 
of HVCs. The absorption line sightlines were not chosen to lie near HI complexes 
(Sembach et al., 2003; Shull et al., 2009), but in these maps it can be seen that 
the majority of the absorbers are in the spatial and kinematic vicinity of the HI 
complexes. In particular, an extension of the Magellanic Stream is evident with 



absorbers at high negative velocities (see also Nidever et al., 2010 Lockman et al., 
2002 Gibson et al., 2000). There are also extensions of Complexes C, A, and 
M and the W complexes, and many of the distance constraints for the ionized 
HVCs are consistent with those of the closest HI HVCs in position- velocity space. 
This map clearly shows that large scale multi-phase flows are active in the halo. 
There are some cases where the absorbers do not show a position-velocity link to 
complexes detected in HI. Most of these are at positive velocities and b > 45°. 
These largely ionized structures could be longer extensions of the W and LA 
complexes, gas along local galaxy structures (i.e., toward the Local Supercluster 
in the north or toward the Sculptor Group for the positive velocity absorber near 



the South Galactic Pole), or outflows from the disk (e.g., Keeney et al., 2006 
Zech et al 



2008) 



2.3 Hot Halo Gas 



The hot gas (~ 10 6 K) in a spiral galaxy halo traces feedback mechanisms ( g4.2| ) , 
and potentially shock- heated IGM (§ |4.1[ ). The extent and mass of this medium 
is of great interest to cosmological models of galaxy formation, as they predict 
it extends out to the virial radius and potentially hosts a large percentage of 



the galaxy's baryons (White & Rees, 1978 


Mailer k Bullock, 2004 


Fukugita & 


Peebles, 2006 


Kaufmann et al., 2006 


Sommer-Larsen, 2006 


Crain et al., 2010 


)■ 



A hot Galactic halo and/or disk- halo medium is inferred to exist from x-ray ob- 



servations in emission (Kerp et al., 1999 



(O VII and O VIII; Williams et al., 2005 



Wang Sz McCray, 1993) and absorption 



measures of pulsars at a range of distances (Taylor et al., 1993 Gaensler et al 



Wang et al., 2005), and from dispersion 



2008 ) . Most of the hot gas detected directly in emission is thought to be within a 



few kpc of the disk (Fang et al., 2006; Yao & Wang, 2007), and this is consistent 



with other galaxies (see £3.3). There is evidence for a lower density medium that 
extends out to at least the Magellanic System (~ 50 — 100 kpc) from several 
indirect methods. 1) The abundant detections of high velocity O VI absorption 
in the Galactic halo is attributed to the interaction of cool halo clouds with the 
hot halo medium (Sembach et al., 2003). 2) Satellite galaxies are stripped of 



their gas out to large radii, most likely as they move through this diffuse hot halo 
medium (IGrcevich & Putman, 20091 iBlitz & Robishaw, 20001). 3) The spatial and 



kinematic structure of the cold halo clouds indicate they are both supported and 
destroyed by a hot halo medium (see ^5]). 4) There is a low-density bipolar struc- 
ture of unknown distance detected in x-ray ( |Bland-Hawthorn Ik Cohen, 2003 ) 
and gamma-ray maps (Su et al., 2010) that is consistent with hot gas being sent 
into the extended halo. 

Mass estimates for the hot halo medium are < 10 10 M within 100 kpc of the 
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disk (e.g.,|Collins et al., 2005| | Anderson fc Bregman, 2010]|Yao et al., 2008[), and 



the covering fraction is estimated to be > 60% (Fang et al., 2006 Sembach et al., 



2003 ) . As discussed further in Q the amount of hot gas in the halo is below that 



predicted by cosmological simulations (sec also Bregman, 2007), but its extent 
makes it relevant to studies of large scale Galactic winds and/or shock- heated 
IGM. The precise physical properties of the extended hot halo medium remain 
uncertain, though densities are generally determined to be between 10 -5 and 10 -4 
cm~ 2 at 50-100 kpc with the indirect methods mentioned above (Putman et al., 



2011b| |Hsu et al., 2011] IStanimirovic et al., 20061 IGrcevich & Putman, 2009 



Sembach et al., 2003). The temperature is usually assumed to be 1 — 2 x 10 



K for these halo density estimates, consistent with the x-ray observations and 
with expectations for gas in equilibrium with a Milky Way mass halo. This 
temperature also happens to be consistent with the velocity dispersion of the HI 
HVCs (~100-120 km s _1 ). There is a great deal of additional work required to 
understand this dimcult-to-detect, low density, hot halo medium. In particular, 
the distribution of this medium and its density at the virial radius need to be 
measured in order to obtain a more accurate estimate of its total mass. 



2.4 Dust, Metals and Molecules 



In this section we discuss the observations of dust and metals in halo gas and 
the limited detections of molecules. Halo dust and metals are direct tracers of 
the ejection of material from galaxies as they result from stellar evolution and 
supernovae. Small dust grains are the dominant heating mechanism of the cold 
neutral medium and metals the dominant coolant, so their abundances dictate 
the temperature structure of neutral halo clouds (Wolfire et al., 1995). Dust is 
destroyed by grain-grain collisions and by particle-grain collisions (sputtering) via 
shocks (Jones et al., 1996 Draine &; Salpeter, 1979a|b|), and therefore halo dust 



properties can provide information on the dynamical history of the gas such as the 
launching mechanism. Understanding halo dust is also important for correcting 
the observed colors of extragalactic objects. The standard Galactic extinction 
map was made using a single grain size and composition for all Galactic gas, and 
a single temperature along a given line of sight ( |Schlegel et al., 1998[ ) . Metallicity 
measurements can further help to determine the origin of the gas (see Q , while 
molecular gas may point to future star formation in the halo. 

Dust can be detected in the Galactic halo through far-infrared (IR) emission, 
depletion of refractory elements, and the reddening of light from distant objects. 
While far-IR emission provides a clean indication of dust in the disk, it is diffi- 
cult to detect in the halo given the lower HI column densities and metallicities 
of halo clouds compared to the disk. HVCs may have less of their metals in dust 
grains (Wakker et al., 2001) and have cooler (and dimmer) dust. No definitive 
detections of HVC dust have been made with far-IR emission, although there 



are tentative detections of Complexes M and C (|Peek et al., 2009 


laboration et al., 2011 


Miville-Deschenes et al., 2005 


). Measuring 



of refractory elements onto grains can in principle indirectly show the presence 
of dust, but it is difficult given the unknown ionization states of the gas across 
species, the likely variation in inherent abundance patterns for clouds, and even 
the significant disagreement in the baseline non-depleted solar abundance ratios 
(Asplund et al., 2009). Nevertheless, there is no sign of the depletion of S and 



Si onto grains for Complex C (Collins et al., 2007 Tripp et al., 2003) or ionized 
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HVCs ( Richter et al. , 2009 ) , while evidence for dust depletion has been found for 



clouds associated with the Magellanic System (Lu et al., 1998, Fox et al., 2010) 



Reddening of background sources is independent of temperature and ionization 
state and may be used in the future to detect or rule out dust in Galactic halo 
gas at interesting levels. 

Metallicity measurements trace the star-formation history of halo gas more 
directly than the dust. Metallicities have been measured primarily only for the 
large Complex C and Magellanic System complexes (MS and LA; see Figures [I] 
- [3]) , although additional estimates will be obtained with future Cosmic Origins 
Spectrograph (COS) observations (Wakker, 2009). Metallicities for Complex C 
are 10-30% solar, and there is some evidence for mixing with Galactic material 



at low latitudes (Shull et al., 2011; Collins et al., 2007). The nearby Complex A 



has a similar measurement (Wakker et al., 2001). The deuterium abundance was 



also measured for one sightline in Complex C and is consistent with cosmological 
expec tations (N(D I) = 2 x 10 15 cm" 2 ; D/H = 2.2 x 10~ 5 ; |Sembach et al., 
2004). The Magellanic Stream and Leading Arm metallicity determinations are 
also ~ 10 - 30% solar (|Fox et al , 20101 IGibson et al., 20001 |Lu et al., 1998 



Sembach et al., 2001). This metallicity, and the relative abundances, have been 



used to support an SMC origin for MS and LA; however, the metal content of 
these features depends on when, and from what region, the gas was stripped from 
the Magellanic Clouds. 

Molecular gas has only been detected along a few sightlines through HVCs. 
H2 was detected in the Leading Arm and Magellanic Stream (~ 10 16-17 cm 



-2. 



Sembach et al., 2001, Richter et al., 2001). It was not detected toward a sightline 



through Complex C, but this may be due to a low HI column along the sightline 
(Richter et al., 2003). CO and H2 was detected in the bridge of gas connecting 



the LMC and SMC, which is consistent with the star formation in the Magellanic 



Bridge (Muller et al., 2003 Lehner, 2002). 



2.5 The Disk-Halo Interface 

Observations of gas at the disk-halo interface help to determine the relationship 
between halo gas that comes from the disk and halo gas that fuels the disk. For the 
Milky Way, the gas at the interface between the disk and halo has traditionally 
been studied separately from the halo gas. The notable exception is the hot 
gas because the velocity resolution of these data do not allow the two to be 
separated. Defining individual objects at the disk-halo interface is difficult at 
any wavelength. Objects are defined with HI data with the large intermediate 
velocity clouds (IVCs) and a population of small discrete clouds, but at other 
wavelengths the properties of the entire gas layer are usually presented. 

For the cold gas, the thickness of the HI layer in the inner Galaxy is 220 pc, 
but approaches < 100 pc in the nuclear regions and flares dramatically at R 
> Rq (Dickey &: Lockman, 1990 see Figure k3J). In some regions, there are kpc 



extensions from the disk due to the presence of IVCs or the Galactic warp ( Levine 



et al., 2006; Kalberla & Dedes, 2008). There is observational evidence that the 
HI temperature increases with increasing z-height. This is also seen in the warm 
and warm- hot gas (see below). Determining the variation in the movement of the 
gas with z-height is difficult given our location. The model of Levine et al. (2008) ) 
has a vertical falloff of the rotation curve of -22 ± 6 km s _1 kpc _i near the disk, 
while Marasco & Fraternali (2011) find lower values of -15 ± 4 km s _1 kpc -1 . 
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Both the temperature and kinematic gradient of the disk-halo gas is depicted in 
Figure [5j 

The large HI IVCs are not easy to isolate as individual objects due to their 
connection to Galactic HI. They are found to reach z-heights of 1-2 kpc and 
their HI mass estimates are 0.5 — 8 x 10 5 M . IVCs have near-solar metallicities 



(Wakker et al., 2001, 2008), and dust is readily detected through IR emission. 



IVCs tend to have somewhat smaller, hotter dust grains than the disk, consistent 



with dust shattering through shocks (Heiles et al., 1988, Peek et al., 2009 Planck 



Collaboration et al., 2011). Small amounts of molecular gas have been detected 



in numerous IVCs with H2 column densities of 10 14 17 cm 2 , consistent with the 



presence of dust (Richter et al., 2003). Deuterium has also been detected (Savage 



et al., 2007). 



Recently a large number of small, discrete HI clouds have been identified as 
a distinct disk-halo population and modeled as clouds rotating with the Galaxy 



at z-heights of 1-2 kpc (|Begum et al. 2010} |Saul et al., 2012| |Ford et al., 2010 
Lockman, 2002} |Kalberla fe Kerp, 2009| ). The GALFA-HI survey has classified the 



largest population of these clouds (~ 1250), and can split the discrete HI clouds 
at intermediate velocities (20-90 km s _1 ) into a cold (T < 5000 K) population 
that appears to be relatively uniformly distributed and related to the disk, and 
a warm population that has a bias toward negative velocities and is distributed 
similarly to the large IVCs (Saul et al., 2012). The HI mass of the clouds depends 
on their distance and range from ~ 1 Mq for clouds at high Galactic latitudes 
at a distance of 1-2 kpc (Saul et al., 2012), to ~ 700 M Q for clouds toward 



the inner galaxy with distances obtained with the tangent point method (Ford 



et al., 2010). The metallicity and dust content of these clouds has not yet been 



measured. 

The warm gas detected in Ha emission at the disk-halo interface is referred 



to as the Reynolds or warm ionized medium (WIM) layer (Haffner et al., 2003 
Reynolds, 1993). This is a layer of warm gas extending 

0.1 cm 



with a volume averaged density of 0.01 



at 1-1.5 kpc (Reynolds, 1991; Gaensler et al., 2008 Haffner et al, 2009). 



2 kpc above the disk 
3 and a filling factor of > 30% 

As 



illustrated in Figure |5j at lower z-heights more of the volume is filled in with HI 



and at higher z-heights the hot gas fills more of the volume (Savage & Wakker. 



2009; Gaensler et al., 2008). The scale height of gas at multiple temperatures, 



as well as its relation to the Galactic disk, has been investigated via absorption 
line studies (Savage h Wakker, 2009 Bowen et al., 2008). The cool and warm 
components have scale heights calculated from the absorption line results that 
are consistent with that obtained from the emission line observations. The warm- 
hot gas, that is often referred to as transition temperature gas, is observed with 
Si IV, C IV and O VI absorption lines and has a larger scale height (~3-4.5 kpc) 
than the cold and warm gas. This warm-hot halo gas is found toward the vast 
majority of the sightlines in the Galactic disk. The disk-halo interface region is 



discussed further in £6.2 



3 Halos of Other Spiral Galaxies 

Halo gas observations for other spiral galaxies provide insight into the unique- 
ness and location of the Milky Way halo gas and the likely origin of the gas. 
Extragalactic observations are generally not as deep as Milky Way observations 
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in the case of emission line observations, and limited to a single sightline for 
absorption line observations. For emission line observations, the break between 
what is disk gas and what is halo gas depends on the spatial and kinematic model 
of the galaxy adopted. Because the halo gas usually connects to the disk gas in 
extragalactic observations, the break is often somewhat arbitrary and varies in 
the literature. For absorption line observations, we are limited to the projected 
distance and velocity of the absorber as a means to associate the gas with the 
halo of a galaxy. 

This section provides a brief overview of observations of halo gas in other 
galaxies with the same sub-divisions as ^2} The exception in this section is 
that there is no attempt to discuss disk-halo gas separately, since the physical 
resolution of extragalactic observations is generally poorer. The spiral galaxies 
noted in this section are just a few examples, and though they are generally Milky 
Way- like, they do span a range of total mass and SFR. Differences in intrinsic 
galaxy properties can lead to different expectations for the amount of halo gas 
(see @. 



3.1 Neutral Hydrogen Halo Gas 



HI halo clouds have now been detected beyond the disks of numerous spiral 
galaxies. The lowest HI mass clouds observed around other spirals is ~ 10 5 



Mq (Thilker et al., 2004), but smaller clouds are likely to be detected with deeper 
and higher resolution observations. The detected clouds are generally within 
~ 10 kpc of the disk, and deep observations do not show evidence for HI clouds 
at distances > 50 — 80 kpc f rom the disks (Westmeier et al., 2007; Pisano et al., 

w 



2007). 



Sancisi et al. (2008) find that extraplanar features with masses > 
25% of the field galaxies. 



Mq exist in the halos of ~ 25% of the field galaxies. 10 M Q is a very large 
HI feature (on the order of the Magellanic Stream), and if the criteria for a HI 
feature to be distinct is relaxed to include warps and/or kinematic evidence for 
deviations from disk rotation, the detected fraction increases to ~50% (Sancisi 



et al., 2008 Haynes et al., 1998). Thus far, the majority of the spiral galaxies 



observed to deep enough levels show a disk-halo component that lags the rotation 



of the disk (often referred to as an 'HI beard'; Sancisi et al., 2001) with a typical 
gradient of magnitude 15-30 km s _1 kpc -1 (Fraternali et al. 



2002 



Oosterloo 



et al., 2007 Heald et al., 2011). These deep observations also often find very 



limited amounts of distinct halo gas (Irwin et al., 2009; Zschaechner et al., 2011) 



Many of the halo HI features that can clearly be disentangled from the disk 
of the spiral are potentially linked to satellite accretion or interaction with a 



companion 


- e.g., M31 f 


Westmeier et al., 2005a 


Ibata et al., 200'i 


7 ), NGC 891 


(Oosterloo et al., 2007 Mouhcine et al., 2010), M33 (Putman et al, 2009c Grossi 


et al., 2008 


), NGC 5055 


(Battaglia et al., 2006 1 


Vlartmez-Delgado et al., 2010), 


NGC 253 ( 


Boomsma et al., 2005 


Beck et al., 1982), NGC 2442 (] 


ilyder et al., 


2001); and see Hibbard et al. (2001 


) and Sancisi et al. (2008 ) for other examples. 



for the Milky Way with only the Magellanic System HVCs being distinct if it 
were observed externally. Future deep stellar observations may link additional 
HI halo features to satellite accretion; however, this is not possible when the HI 
feature overlaps with the disk in projection - e.g., NGC 2403 (Fraternali et al., 



2002, 2001); NGC 4395 (Heald et al., 2007). 



The high resolution observations of galaxies designed to detect halo clouds 
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often have a HI column density sensitivity of ~ 10 cm . This is only the 
typical peak column density of a Milky Way HVC and therefore may be missing 
a significant amount of halo gas. The spiral galaxy that has been mapped to the 
deepest levels (N#j < 10 18 cm -2 ) at a resolution of a few kpc, and is most similar 



to the Milky Way, is M31 ( |Thilker et al., 2004[ fWestmeier et al., 2005a| ). NGC 891 
has also been mapped to 10 i8 cm -2 levels at this resolution, but this galaxy has 



a much higher SFR than the Milky Way (Oosterloo et al., 2007). Both of these 
galaxies are shown in the top panels of Figure [6} The majority of their halo gas 
is connected to the HI disk (like most other galaxies), and/or could be associated 
with a satellite interaction. The calculated HI mass in M31's halo is the same as 
that calculated for the Milky Way without the Magellanic System (3 — 4 x 10 7 
M ; see £2.1.1), and presumably these clouds also have a substantial ionized 



component similar to the Milky Way's HVCs (Fittingoff et al., 2009). NGC 891 



shows a larger amount of halo and disk-halo gas, which is consistent with its 



higher star formation rate (Oosterloo et al., 2007). A diffuse halo medium can be 



inferred to exist from the head-tail nature of several of the M31 HVCs (Westmeier 



et al., 2005a). The M31 HVCs also have similar linewidths to Galactic clouds. 



3.2 Warm Ionized and Warm-Hot Halo Gas 

Similar to the Milky Way, the warm gas in the halos of other spiral galaxies is 
detected in emission and absorption. The emission line observations primarily 
detect gas close to the disk of the spiral (within ~5 kpc; Hoopes et al., 1999 



Collins Sz Rand, 2001 ), and show that the extent of the extraplanar Ha emission 



scales with the galaxy's SFR (e.g., Rand, 1996 Rossa & Dettmar, 2003). Invest- 



igations of the motion of the ionized gas find a similar rotation lag to the HI with 



increasing z-height ( |Rand, 2000} |Heald et al., 2006| ). |Heald et al. (2007[ ) find a 
typical vertical gradient of magnitude 15-25 km s _1 per scale height. The gas 



increases in temperature with z-height (Collins Sz Rand, 2001 Rand et al., 2011 ) 



as also found for the Milky Way. See Haffner et al. (2009 ) for a further review of 



the warm ionized gas of spiral galaxies detected in emission. 

There are abundant absorption line results that are presumed to trace the 
warm gas in galaxy halos based on their observed kinematic and spatial proximity. 
Since low column density absorbers (N#/ ^10 15 cm -2 ) are also detected along 
large scale cosmic filaments, the association with a galaxy halo is not certain 
(IBowen et al., 19961 IPutman et al., 20061 IPenton et al., 20021). Within ~ 300 kpc 



of a galaxy, Lya absorber results sh ow a covering fraction close to 100% and total 
gas mass estimates are 10 9-10 M Q (Bowen et al., 2002 Wakker Sz Savage, 2009 



Prochaska et al., 2011; Chen Sz Prochaska, 2000). For the higher column density 



Lya absorbers (and Mg II absorbers that are generally thought to trace gas with 
Nhi > 10 16 cm -2 ), a relationship to the halos of galaxies is relatively clear as 



the detections are almost solely in the vicinity of galaxies (Bowen et al., 2002 
Penton et al., 2002| Chen et al., 2010a). There are a limited number of these 
high column density absorbers at z = around spiral galaxies. Studies of Mg II 
absorbers at z ~ 0.1 — 1 suggest the covering fraction is > 50% within 100 kpc 



of the galaxy's center (Chen et al., 2010a Kacprzak et al., 2008); although the 



covering fraction is expected to be higher at higher redshifts (Fernandez et al., 



2012) 



Studies of the warm-hot component (~ 10 5 6 K) are largely limited to O VI 
absorbers at z = (iProchaska et al., 20111 IWakker Sz Savage, 2009), but will 
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soon be extended to include C IV (Green et al., 2012). Broad Lya absorbers 



have also been claimed to be warm- hot gas (Bowen et al., 2002; Savage et al., 



2011 Danforth et al., 2010). The low redshift O VI absorbers have significant 



covering fractions of ~ 70 — 80% within 200-400 kpc of galaxies (Prochaska et al., 
2011 Wakker &: Savage, 2009| ). This seems to be consistent with the extensive 
work at higher redshifts that measure covering fractions close to unity for C IV 
{z ~ 0.5) and O VI (z > 0.1) within 100 kpc of star forming galaxies (Chen 



et al., 2001; Tumlinson et al., 2011 Chen & Mulchaey, 2009). Star forming 



galaxies are more likely to have both warm and warm-hot gas in their halos 



( Tumlinson et al, 2011] |Chen et al., 2010b[ ). |Tumlinson et al. (201 i[ ) estimate 
a mass of circumgalactic gas of 2xl0 9 M (assuming solar metallicity) out to 
150 kpc from the abundant O VI detections near galaxies at redshifts z = 0.1- 
0.36. This is similar to the mass derived from the Lya absorber detections in 
galaxy halos noted above, but approximately an order of magnitude higher than 
the mass estimates for ionized gas in the Milky Way halo (£ |2.2[ ), This may be 
due to several reasons: 1) extragalactic absorption line observations may also 
capture gas in cosmic filaments near the galaxy, but beyond the virial radius, 2) 
the assumptions used to calculate total gas mass from extragalactic absorption 
line results lead to an over-estimate, 3) a range of galaxy types are used for 
extragalactic measurements, or 4) the Milky Way warm and warm- hot halo gas 
is less massive or has not yet been fully observed. 



3.3 Hot Halo Gas 

The x-ray detected hot gas in the halos of Milky Way-like spiral galaxies always 



extends from the galaxy's disk and is within a 10 kpc radius ( 


Li et al., 2008 


2006, Wang et al., 2003; Immler et al., 2003 


Otte et al., 2003 


Wang et al., 2001 



Bregman Sz Houck, 1997). The scale heights are calculated to be 1-2 kpc for 



many spiral galaxies, and the x-ray features are similar to the extraplanar fea- 



tures detected in Ha emission (Strickland et al., 2004; Bregman & Houck, 1997 



Tyler et al., 2004). There are numerous limits on the amount of hot halo gas 



at large radii from non-detections (e.g. Yao et al., 2010 Anderson Sz Bregman, 



2010, Rasmussen et al., 2009). As with the Milky Way, most researchers ar- 
gue that the hot halo medium cannot account for a large percentage of a spiral 
galaxy's baryons. If galaxies substantially more massive than the Milky Way are 
considered, the hot halo gas extends to larger radii (e.g., ~50 kpc for NGC 1961 
with M tnt > 10 13 M ; | Anderson & Bre gman, 201lb. Unlike the Milky Way, there 



is limited indirect evidence for an extended, diffuse, hot halo medium around 



other spiral galaxies (although see {3.2 for the warm-hot gas). One exception is 
M31, where the head-tail HI clouds out to ~ 30 kpc from the disk are likely to 
originate from the movement of cold clouds through a diffuse, hot halo medium 



(Westmeier et al., 2005a). 



3.4 Dust, Metals and Molecules 

Direct detections of dust are primarily at the disk-halo interface of spiral galaxies 
(e.g., |Keppel et al., 1991[ ). A survey of nearby edge-on galaxies found dense 
dust {Ay ~ 1) at least 2 kpc above the plane, and the dust is correlated with 



the presence of extraplanar diffuse ionized gas (DIG) (Howk Sz Savage, 1999). 
Menard et al. (2010| ) have made a convincing detection of dust in extended galaxy 
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halos using the reddening of SDSS background quasars. They found substantial 
reddening 20 kpc from the galaxies, with detectable reddening all the way out 
to 10 Mpc (see also Zaritsky, 1994). Under an assumption of SMC-like dust, 
they claim that galaxies have similar amounts of dust in their halos and disks, 
consistent with observations of reddening in Mg II absorbers at intermediate 
redshift flMenard & Fukugita, 2012] ). 

Metals are also present in spiral galaxy halos, as directly evident from their 
detection with absorption lines ({ 3.2). Though ionization conditions are difficult 
to assess, Tumlinson et al. (2011) estimate there is ~ 10 7 Mq in oxygen within 
150 kpc of star forming galaxies at z = 0.1 — 0.36 using O VI measurements. 
The abundant detections indicate enriched gas has been fed into halos for some 
time. Sightlines for absorption line experiments that pass through extragalactic 
HI halo features are rare, but in the cases where star formation has occurred in 
the gas, metallicities can be measured with the HII regions. Most of this type of 
halo gas is related to an interaction and the metallicity of the halo HI features is 
similar to the outer region of the galaxy (Werk et al., 2011). These interacting 
systems with halo star formation are also the only ones for which molecular gas 
has been detected or inferred to exist (e.g. Lee et al., 2002). 



4 Origins 



Spiral galaxy halo gas is likely to originate from multiple sources, including from 
the intergalactic medium, satellite galaxies, and the disk of the galaxy itself. This 
section discusses these origins in the context of the observational results presented 
in ^2] and §[3j Viable origin scenarios should be consistent with the following key 
observational features. 

• Large scale flows of multiphase g prevalent in galaxy halos. This is 
especially evident from the spatial and kinematic link between the cold, 
warm, and warm-hot gas in the Milky Way halo (see Figures [l] and [2| . 

• The densest halo gas is closest to the disk in all phases. Exceptions (and 
the largest cold halo features) are largely linked to satellite accretion (see 
Figures [3] and [6]). 

• Gas at the disk-halo interface shows a clear temperature gradient (cold to 
hot), and kinematic gradient (decreasing rotation) with increasing z- height 
(see Figure [5j). This gas is also more enriched than gas further out in the 
halo and more extensive in galaxies with higher SFRs. 

• Diffuse gas extends to radii > 100 kpc in spiral galaxy halos. This is 
evident for the Milky Way from the properties of clouds associated with 
the Magellanic System and stripped dwarf galaxies, and for other galaxies 
from absorption line results that show large covering fractions. 

• For the Milky Way, the mass in warm and cold halo gas is < 10 9 M , and, 
though largely model-dependent, the hot gas mass is < 10 10 M . These 
measurements are consistent with the HI and x-ray observations for other 
low redshift spiral galaxies of similar mass. 

• The halo gas appears to be bound to the dark matter halos of the spiral 
galaxies. This gas is also found to be enriched with metals at some level 
(i.e., non-primordial). 



17 



Throughout this section we refer to and make connections to a cosmological 
grid simulation of a Milky Way-mass galaxy at z = (referred to as MW simula- 
tion; Joung et al., 2012a Fernandez et al., 2012). We use this simulation not as 
the final word on the state of the Milky Way halo, but rather as a guide to our 
understanding of the important physical processes at play. The simulation, per- 



formed with an adaptive mesh refinement (AMR) code, Enzo (Bryan & Norman, 



2000 O'Shea et al., 2004), has sufficiently high spatial (136-272 pc comoving or 



better at all times) and mass (~ 10 5 M ) resolutions to allow us to study and 
track the spatial and kinematic distribution of the multiphase gas in the halo in 
detail. The simulation includes metallicity-dependent cooling, metagalactic UV 
background, shielding of UV radiation by neutral hydrogen, photoelectric heating 
for the diffuse gas, star formation and stellar feedback (but not AGN feedback). 
Figure [7] displays a snapshot of the MW simulation in HI at z = 0.3 and zooms 
in on halo gas origins relevant to each of the below sections. Figure [8] shows 
a cartoon based on the simulation that is split to show the main accretion and 
feedback mechanisms distinctly. It also depicts the density-temperature diagram 
of simulation cells at z=0. 



4.1 Accretion from the IGM 



Under ACDM, the large majority of the baryons enter the virial radius of a 
galactic halo in gaseous form along cosmic filaments. This is in contrast to the 
classical picture that posited that the incoming gas would spherically collapse 
and be shock-heated close to the virial temperature of pressure-supported gas 
within a galaxy's halo (IWhite & Rees, 19781 IWhite & Frenk, 19911). Analytic 



arguments (Birnboim & Dekel, 2003 Dekel <fc Birnboim, 2006 Binney, 1977) and 



SPH simulations have shown that some of the gas will never be shock-heated 



to the virial temperature of the halo and be accreted via a "cold mode" (Keres 



et al., 2005; Brooks et al., 2009). This mode is found to dominate in galaxies at 
high redshift (z ^ 2), as well as in present day low- mass galaxies (Mh a io ~ 10 12 
M©). We note that this mode is referred to as cold, but in reality the gas is 
at temperatures < 10 5 ' 5 K, so much of it is not cold under our observational 
definition in ^2| 

The distinction between hot and cold mode accretion is based on the maximum 
temperature that a given gas particle attains along its trajectory as it enters the 
virial radius. At z ~ 0, for galaxies with masses similar to the Milky Way, a large 
fraction of the accreted gas is hot mode and does not connect directly to the 
disk. The hot mode gas is heated either by one strong shock close to the virial 
radius or via a series of weak shocks as the kinetic energy of the outer portions 
of the filaments is converted to thermal energy due to compression from the 
existing halo gas (Keres et al., 2005, Joung et al., 2012a). In the MW simulation, 



~ 70% of the mass inflow at the virial radius is consistent with coming in along 
filaments. The cold mode of the accretion at z ~0 is found in the core regions of 
the compressed filaments at temperatures ^ 10 5,5 K. This gas stays cold because 
it is shielded from interacting with the existing halo medium by the hot mode 
gas that surrounds it in the filamentary flow (Joung et al., 2012a). The accretion 



of the hot and cold mode gas along filamentary flows is depicted in the top panel 
of Figure [8] and box 2 of Figure [7j 

The properties of the filamentary flows can be examined further in the MW 
simulation. In the density-temperature diagram of Figure |8l the filamentary 
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flows generally have intermediate specific entropy between that of the disk gas 
and the outflows. The flows come in at temperatures of ~ ]_q4.5— 5.5 an( j 
outer portions are heated toward ~ 10 6 K and the inner portions cool toward 
~ 10 4 K. The gas at < 10 5 ' 5 K has a relatively high neutral fraction (generally 
fneut > 10~ 4 ), and extends from radii of approximately 25 - 200 kpc. Some of the 
cold cores of the filaments have HI densities comparable to those of gas stripped 
from simulated satellites, but they generally have lower metallicities (Z/Zq < 0.2; 
see Figure [7]). Details on the HI properties of the halo gas in the MW simulation 
are outlined in Fernandez et al. (2012). 

In terms of observations, the filamentary flows could be represented by the 
diffuse warm and warm-hot gas detected with absorption lines near the velocity 



of the galaxy (£3.2); the majority of these filaments are not detectable in HI with 
our current sensitivities until close to the disk (see below). The covering fraction 
of the observed warm halo gas approaches 100%, and this is consistent with the 
covering fraction in the M W simulation when low e nough column density gas is 
selected (< 5 x 10 13 cm~ 2 ; Fernandez et al., 2012). The gas that is heated to 



temperatures > 10 5,5 K largely becomes part of the general hot halo medium, and 
may be related to the fact that we s ee ev idence for a diffuse hot halo medium 
out to large radii (however see also ^4.2). The MW simulation finds 2 x 10 10 
Mq in hot gas within the virial radius and this is consistent with the simulation 
of Sommer-Larsen (2006| ) . This is greater than the amount found by observations 
thus far. 

For the gas in the filamentary flows to feed the star-forming disk, the gas 



needs to cool further; the gas potentially forms clouds within the halo (Mailer & 
Bullock, 2004[|Connors et al., 2006a] |Kaufmann et al., 2006] [2009} [LirTfc Murray, 
2000). The cold mode gas contains density enhancements that can cool as long as 



the cooling time is shorter than the disruption time ( Joung et al., 2012b; Binney 



et al., 2009). Linear fluctuations were not able to cool in these detailed studies 



and cooling without significant density enhancements is also not seen in the MW 
simulation (Figures [7] & [8]) . The gas that does cool may not survive the trip to 
the disk without disrupting further (see ^5]). It remains somewhat unclear how 
halo gas accreted in hot mode can feed the disk. Possibly closer to the disk, where 
it is denser, the hot gas may mix with feedback material and become part of the 
disk-halo interface (£ |6.2[ ). In any case, in the MW simulation the gas is only 
observable as HI clouds within ~ 50 kpc, and this is consistent with observations 



and previous simulations (Sommer-Larsen, 2006; Kaufmann et al., 2006 
Hernquist, 2009). The mass in observed warm and cold clouds is 10 s_ 
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both the observations and simulations when the same type of gas is selected ( Peek 



et al., 2008; Fernandez et al., 2012). 



4.2 Feedback 

After baryons collect at the centers of galaxies and form stars, the stars and 
central black holes act back on their environments and can affect the gaseous 
galactic halos. These "feedback" processes include radiation from young stars 



that heat and ionize the halo gas (Bland-Hawthorn & Maloney, 1999, 2002, Bland- 



Hawthorn & Putman, 2001), mechanical energy from supernova explosions (e.g., 
Mac Low & Klessen, 20041 IVeilleux et al., 20051) and AGNs (lAntonuccio-Delogu k 



Silk, 2010 Di Matteo et al., 2005, 2008), as well as momentum-driven winds due 



to radiation pressure from central sources on dust grains (Murray et al., 2005) 
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and a galactic fountain as enriched hot gas rises up from the disk (e.g., Shapiro & 



Field, 1976; Joung & Mac Low, 2006). We will discuss each of these mechanisms 



in turn. Feedback is required in cosmological simulations to reproduce the funda- 
mentals of observed galaxies; however, despite its importance, it is described by 
phenomenological prescriptions partially due to our lack of understanding of the 
dominant process(es) (Oppenheimer & Dave, 2006 Springel & Hernquist, 2003 



Stinson et al., 20061 Dalla Vecchia Sz Schaye, 2008). Halo gas is in many ways a 



more direct consequence of disk feedback than the stellar observables of galaxies 
and is therefore an important tool for understanding these processes. 

Radiation escaping from a galaxy's disk ionizes and heats halo gas. Ha obser- 
vations of cold Galactic halo clouds have been used to estimate that 6% of the 
Milky Way's ionizing photons escape normal to the disk (f esc = 1-2% when av- 
eraged over the solid angle; Bland-Hawthorn &: Maloney, 1999] , and see { 2.2 and 



£2.1.1). For other spiral galaxies, there are mainly upper limits (e.g., Tumlinson 



et al., 1999 Deharveng et al., 2001), and the measurements are consistent with 



escape fractions of less than a few percent as long as starbursts are excluded. 
As expected, if f esc is estimated with a galaxy's diffuse ionized gas it is highly 



dependent on the dumpiness of the galaxy's ISM (Clarke & Oey, 2002; Zurita 



et al., 2002 Fernandez &: Shull, 201 lj ). Radiation escaping from galaxy disks is 



only starting to be included in cosmological galaxy formation simulations. Con- 
straints on the escape fraction impact both halo gas and the ionization of the 



IGM as a whole, in particular at higher redshift (e.g., Ciardi et al., 2003 Gnedin, 



2000). 



Mechanical feedback mechanisms, such as SN-driven or past AGN-driven out- 
flows, are the most relevant in the context of the origin of a spiral galaxy's halo 
gas. For the Milky Way, a diffuse hot halo medium is evident out to the distance 
of the Magellanic System (~ 50 — 100 kpc), and this indicates either relatively 
strong Galactic winds and/or gas left from hot mode accretion is present. The 
detected x-ray and gamma-ray 'bubbles' are also consistent with the Milky Way 
having feedback from a central engine in the past 10-15 Myr (Su et al., 2010). 
For other galaxies, outflows are consistent with the enriched warm- hot gas at 
distances on the order of the galaxy's virial radius (§ |3.2[ ). Finally, the baryon 
fraction within the virial radius of Milky Way-like spirals is substantially smal- 
ler than the cosmic mean (~ 0.4 for the Milky Way), and this is attributed to 
powerful outflow mechanisms (Bregman, 2007 |McGaugh et al., 2010 ), as is the 
ability to reproduce the basic properties of a galaxy in cosmological simulations 

2MT| . 



(Benson et al., 2003; McCarthy et al., 2010 



Some of the material from large scale galactic winds is likely to escape the 
galaxy's halo (potentially explaining the lack of baryons), and a fraction is also 
thought to feed the hot halo. In Figure [8] we show a cartoon based on the MW 
simulation that depicts outflow due to SN-driven winds. This wind extends from 
velocities of 200 km s -1 to over 1000 km s _1 . The fastest gas escapes, while the 
lower velocity gas feeds the hot halo high entropy, enriched material. Observa- 
tionally, the lower velocity material may be some of the low density absorption 
line systems in the Milky Way and other galaxy halos (^2.2| &; {3.2); in partic- 



ular, the gas at high positive velocity in the general direction of the Galactic 
Poles shown in Figure [2] The SN-wind gas at higher velocities is too hot and low 
density to be observable. Figure [8] also shows that the MW simulation confirms 
previous simulation expectations that supernova feedback can be thought of as 
injecting high entropy gas into the system. In particular, once launched, the 
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specific entropy of this enriched gas remains unchanged for a large range of the 
galactocentric radius, suggesting adiabatic expansion. 

Momentum driven feedback is also commonly discussed as a method of launch- 
ing material from the disk. This mechanism is likely to be relevant only in galaxies 



with powerful starbursts or AGNs (Oppenheimer <fe Dave, 2006; Murray et al., 



2007) that may have operated in the more violent, formative years of present-day 
spiral galaxies like the Milky Way (z ~ 2 — 3; e.g., Kim et al., 2011 ). In addition, 
this mechanism may have difficulty launching material as, even if one assumes 
the maximal efficiency of converting mass to energy and that all the available 
momentum is used to drive outward motion of the surrounding gas, the outflow 
velocity would be too low to eject gas for reasonable values of the mass loading 
factor. Momentum driven winds have also been invoked as a method to launch 



substantial amounts of dust into galaxy halos (Aguirre et al., 2001); however, 
energy-driven mechanisms may be equally effective through the coupling of dust 
and gas. Dust observations thus far could be consistent with either mechanism. 

A galactic fountain is distinct from the powerful galactic winds discussed above 
in that it is a disk-wide phenomenon in which gas is driven upward by the cumu- 
lative effect of supernova explosions and adiabatic expansion, and subsequently 



falls back down due to radiative cooling (Shapiro & Field, 1976 Bregman, 1980) 



The galactic fountain phenomena is most likely to affect the detailed structure 
and dynamics of gas at the disk-halo interface. Kiloparsec-scale ISM simulations 



in vertically elongated boxes ( 


de Avillez, 2000; 


de Avillez & Breitschwerdt, 2005 


Joung & Mac Low, 2006 


Booth & Theuns, 2007; 


Hill et al, 2012 


) report ejecta 



cooling and condensing out at a few kpc heights having morphology and kin- 
ematics similar to those of observed IVCs (£ |2.5[ ). The near-solar metallicities and 
dustiness of many IVCs are also consistent with a disk origin. The temperature 
gradient of the disk-halo interface gas is expected with hot gas from the SN rising 
to higher z- heights, and the kinematic lag suggests a combination of this type of 
feedback and fueling (see £|6.2[), The large population of small disk- halo clouds 



may also be linked to star formation processes in the disk (Ford et al., 2010). The 
cartoon shown in Figure [5] summarizes some of the observational findings for gas 
at the disk-halo interface and has many similarities to the fountain simulations. 

Feedback material mixes and interacts with the existing halo gas. In box 1 of 
Figure [7j we show an example in the MW simulation of neutral halo clouds that 
have originated from the interaction of a supernova-driven wind and incoming 
cooler flows. The gas has condensed at distances of 50-100 kpc from the disk and 
has metallicities around 0.3 Zq, which is greater than many of the other simulated 
HI halo features, and lower than that of the wind itself. The HI clouds fragment 
further and impact the disk ~300 Myr later than the timestep shown. The higher 
metallicities of some regions in HVC Complex C may be representative of this 
mixing process closer to the disk. HVC interaction with feedback material may 
also be evident in the Ha brightness of some complexes and the mixture of photo- 
and collisional ionization required to model the line ratios measured for halo gas. 



4.3 Satellite Accretion 



Satellites lose their gas as they move through a spiral galaxy's halo (e.g., Putman 



et al., 2003b Grcevich & Putman, 2009). The gas is primarily ram pressure 



stripped from the dwarf in models, although tidal forces also play a role (Mayer 
et al., 20071 120061). ACDM simulations find 40% of L* galaxies host a 0.1 L, 
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satellite within its virial radius and this is consistent with results from SDSS 
(Tollerud et al., 2011 Liu et al., 2011). This percentage drops significantly at 



smaller radii (e.g., 12% within 75 kpc) and also when only blue and star forming 



satellites like the LMC are considered (8.2% within 100 kpc; Tollerud et al., 2011 



James & Ivory, 2011). We know from the Local Group that the vast majority of 
the satellites have < 0.1 L*, however the amount of gas they bring is substantially 
smaller (~ 10 5-8 M each). For instance in the Local Group, the Magellanic 
System contains ~ 30 — 50% of the HI mass of the Milky Way, and M33 will 
eventually provide ~ 25% of M31's HI mass ( |Putman et al., 2003b[ |2009c[ ). The 
lower mass satellites within the virial radius of the Milky Way will only have 
provided ~ 1%, or if we include the large number of yet unobserved satellites 
predicted by ACDM, ~ 10% ( |Grcevich fc Putman, 2009] ). 

The cumulative effect of the accretion of satellites will leave behind a large 
amount of gas that will eventually feed the disk. Since most galaxies are likely 



to be stripped at radii > 20 kpc (Grcevich & Putman, 2009, Mayer et al., 2006 



Nichols &; Bland-Hawthorn, 2011), much of the cold gas will be integrated into 
the diffuse galactic halo before eventually recooling at the disk (see fj5|. This is 
consistent with multiphase flows being found throughout galaxy halos and helps 
to sync the angular momentum distribution of the satellite gas with the spiral's 
halo gas. This integration of satellite material with the halo medium can be 
seen in box 3 of Figure [7] which depicts the accretion of a triple system of dwarf 
galaxies. The stripped gas shows a temperature gradient as it integrates into the 
hot halo material, and the column density and metallicity of the gas also gradual 
decreases (see also Figure [8]) . The cold and warm satellite gas has on average 
higher radial velocities than other features in the halo. This material becomes 
untraceable several hundred Myrs after being stripped from the satellite core, 
but may leave behind overdensities that help seed the cool filaments closer to the 
disk (see ^5] and ( ggfr . 

The system shown in Figure [7] has many similarities to what is observed for 
the Magellanic System of the Milky Way, which has given us a detailed look 
at the relationship between gaseous satellites and a spiral galaxy's halo gas (see 
Figures [T] - [3] and [6]) . The Magellanic System shows a column density gradient 
along its length and multiphase gas surrounding it, particularly at the tail where 
it is being integrated into the hot halo medium. This system can also account 
for the highest negative and positive velocity halo gas in the Milky Way halo. 
The interaction of the satellite gas with an existing halo medium, and the high 
covering fraction of the warm and warm-hot gas does make it relatively clear that 
mechanisms besides satellite accretion are also feeding the halo. In addition, for 
the Milky Way, most of the HVCs (besides those associated with the Magellanic 
System) do not show any clear link to satellites or stellar streams. There is a large 
body of literature describing the Magellanic System and it is thought to be the 
result of a combination of a tidal interaction between the two satellites and their 
subsequent interaction with the Milky Way ( Connors et al., 2006b| iMastropietro 



ct al 


, 2005; 


2008 


)• 



1999 Diaz fc Bekki, 2012[ [Nidever et al., 



In general if the HI gas observed in galaxy halos is directly from a satellite, 
it is likely to have been relatively recently stripped from the satellite. This is 
consistent with the link between large HI features in spiral galaxy halos and a 
relatively recent interaction event (see { 3.1 ). Satellite HI features are also likely to 



have higher metallicities than gas that originates from the IGM (e.g., Fernandez 
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et al., 2012 ), although given the likely mixing of halo gas components, deep maps 



of the stellar components of galaxy halos may remain the best method to identify 



a satellite origin for the gas (e.g., Martmez-Delgado et al., 2010). The Milky Way 
has a large number of HVCs at z-heights < 10 kpc that do not seem to be related 
to stellar halo features. These clouds may be a combination of satellite gas and 
halo gas of other origins that is recooling close to the disk. 



5 Halo Cloud Survival 



Cold halo clouds moving through hot diffuse galaxy halos are unlikely to survive 
for periods greater than a few hundred Myr unless a strong support mechanism 
is invoked (see below). These short lifetimes are a problem unless the clouds are 
continually regenerated on short timescales. In addition, some HVCs that are not 
clearly linked to satellite accretion have velocities that cannot be reached through 



gravitational acceleration within several hundred Myr (Benjamin & Danly, 1997 



Peek et al., 2007). The exact lifetime of a cloud depends on several factors such 



as cloud density, halo density, and velocity, but the total mass of the cloud seems 



to be one of the largest factors increasing their lifetimes (Heitsch & Putman, 



2009; Kwak et al., 2011). This is consistent with simulations that show the halo 
clouds are destroyed primarily via the Kelvin-Helmholtz (KH) instability. The 
characteristic growth time for the KH instability is txH oc X^^Rd/vrei, where 
X is the density contrast between the cloud and the external medium, R c i is a 
characteristic length of the cloud, and v re i is the relative speed between the two. 
The KH instability leads to gas being ablated from the cloud edges and forming a 
tail of lower column density material. Observationally this is evident in the head- 



tail clouds found throughout the Galactic halo ( 


Putman et al., 2011b 


Briins et al., 


2000 


Ben Bekhti et al., 2006 


Westmeier et al., 2005b 


), and potentially some of 



the detailed structure found along the edges of larger clouds ( Peek et al., 2007 
Winkel et al., 2011) as shown in Figure [4j O VI observations are also consistent 



with the interaction of halo clouds with the hot diffuse halo medium (Sembach 



et al., 2003; Kwak et al, 2011 Fox et al, 2006; Collins et al., 2007). 



Possible methods of supporting the clouds to significantly extend their lifetimes 



include the presence of dark matter (Nichols &; Bland-Hawthorn, 2009; Braun & 



Burton, 2000), magnetic fields ( McClure-Griffiths et al., 2010), and/or shielding 



by an extended diffuse gaseous component (e.g., the extended medium surround- 
ing the Magellanic Stream). The presence of dark matter extends a gas cloud's 



lifetime to timescales on the order of a Gyr (Quilis & Moore, 2001); however 
despite the advantage of this preservation there is no strong evidence for dark 
matter in HVCs. Rather, even the small compact clouds are associated with 



large gaseous complexes (Figure \\\ |Saul et al., 2012 Putman et al., 2011b), and 



it would be difficult for gas to survive in such small dark matter halos through 
reionization (e.g., Ricotti et al., 2008|). A magn etic field may help to preserve the 



clouds (Konz et al., 2002 Santillan et al., 2004 Kwak et al., 2011), however most 



of the simulations done thus far are in two dimensions and some simulations 



indicate a strong magnetic field can aid in the disruption (Stone Sz Gardiner, 



2007). 



Dynamical shielding by outer envelope layers may be the most likely mechan- 
ism of extending the lifetimes of HI clouds. As shown in Figures [7] and § the 
MW simulation finds cold neutral clouds are embedded within larger, warm, low 



23 



density structures and Figure [T] shows that this is also observed. In this config- 
uration, the HI clouds experience smaller relative velocities with respect to their 
immediate surroundings, which increases txH- Head-tail structures or features 
indicative of shearing may be less common in the HI component of those clouds 
that are dynamically shielded by an outer, warmer layer. This may explain why 
clouds associated with the Leading Arm of the Magellanic System do not have 
detectable warm gas and also show abundant head-tail structures, while the tail 
of the Magellanic System has abundant warm diffuse gas and limited head-tail 



structures (e.g., Putman et al., 2011b). 

For the Milky Way, since many of the HI halo clouds are at distances of ~ 10 
kpc (or greater in the case of the clouds associated with the Magellanic System), 
and since the individual clouds within the complexes usually have masses < 
10 5 Mq, the clouds are unlikely to make the trip through the Galactic halo 
unless significantly shielded by an outer, warmer component. For other galaxies, 
the detected large HI clouds are also likely to be collections of smaller clouds. The 
disrupted HI clouds will become part of the multi-phase halo medium, and any 
remaining over-densities will sink toward the Galactic disk. Some of the warm 
and warm-hot gas detected with absorption line systems may be the leftover 
over-densities from disrupted HI clouds. Whether these over-densities are able to 
re-cool depends on the magnitude of the over-density, and, if they cool, they may 



disrupt again before impacting the disk ( Joung et al., 2012b Vietri et al., 1997) 



The infall of these clouds can dissipate energy and help to heat the extended hot 
halo medium (Murray & Lin, 2004). According to multi- wavelength observations, 



gaseous halos are complex multiphase places, which is consistent with an ongoing 
cycle of halo cloud destruction and cooling. 



6 Feeding Galaxy Disks 

The amount, metallicity, and ionization state of the halo gas has consequences 
on a galaxy's resulting star formation. Star formation rates of spiral galaxies at 
low redshift are approximately 0.5 - 5 M /yr, and Galactic chemical evolution 
models find that at least half this amount needs to be continuously accreted, low 



metallicity fuel (Chiappini, 2009 Chiappini et al., 2001). Without accretion, the 



SFRs indicate the star formation fuel currently in the disk will last for another 
few Gyrs at most. In this section, potential sources to maintain a galaxy's star 
formation over time are discussed. This includes a calculation of the accretion 
rate for Milky Way halo gas, a discussion of gas at the disk-halo interface, and a 
brief mention of some often neglected disk sources. 



6.1 Halo Sources 

The mass inflow rate to the Galactic disk is a crucial parameter to derive from 
the physical properties of the HVC complexes outlined in £2.1.1 Many studies 



estimate the HVC contribution to the accretion rate using some form of the 
equation, 

where Mj is the mass of a given complex, Vi is the observed radial velocity and 



Di is the distance (Thorn et al., 2008, Wakker et al., 2007). The largest difficulty 



with this method is that the observed radial velocity of a HVC complex is the 
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sum of the velocity of the complex toward the disk (or more generally a poloidal 
component), and the azimuthal speed of the cloud around the disk (toroidal 
component), projected along the line of sight. This observed velocity may be 
very different from the true radial velocity of the gas toward the Galaxy, and the 
cloud may be accelerated in the future by MW gravity and/or decelerated by ram 
pressure. Additionally, the distance to the gas is not the same as the distance 
from the cloud to where it will eventually impact the disk. 

Now that most of the complexes have direct distance constraints, we can con- 
struct a very simple model of the HI accretion process that takes into account the 
bulk of the accretion onto the Galaxy from known HVCs in the next ~ Gyr. We 
only consider the 13 complexes with relatively well-known distances, excluding 
the Magellanic Stream and Leading Arm complexes as they are unlikely to accrete 
within the Gyr. We calculate their mass as outlined in §2.1.1 which includes the 
clouds' HI, He, and ionized gas detected in Ha, and use the angular position of 
each complex's flux weighted center. The 23% of the HVCs (by flux) that do not 
have a direct distance constraint are not included here, but this gas would not 
substantially increase the calculated accretion rates. We ascribe a radial velo- 
city, v rj0 bs, to the complex equal to its flux- weighted radial velocity in GSR. We 
then construct a very simple model for the 3-space velocity of the system of all 
HVC complexes. All complexes are given an identical azimuthal velocity, vg, and 
accretion velocity toward the Galactic center, vq- Then, for a given vg and vg, 
each complex has some minimum additional velocity to make up the difference 
between the observed radial velocity and modeled radial velocity. The best fit 
parameters for the model are those that minimize these differences. We note that 
our results are qualitatively unchanged if the accretion vector is assumed to be 
toward the disk instead of the Galactic Center. 

The complexes are best fit by [vg,vc] = [77, -40] km s . We note that the 
azimuthal velocity, 77 km s _1 in the direction of Galactic rotation, is roughly 
cons istent with the v r gradient of -20 km s _1 kpc -1 above the disk noted in 
£2.5, and that the accretion velocity of -40 km s _1 is roughly consistent with the 
average radial velocity of all HVCs (~ -50 km s _1 ). We show the accretion of 
complexes as a function of time given this best fit velocity in Figure [9| as well as 
a 'maximal' accretion rate using the highest vq allowable with twice the residual 
in the fit, -200 km s" 1 . The fit and maximal accretion rates are 0.08 and 0.4 
M /yr, respectively. The maximal value approaches that required by chemical 
evolution models for the Milky Way (at least 0.45 M /yr; Chiappini, 2009), but 
the fit is far below that required. These values are consistent with accretion rates 



of neutral gas from simulations (Peek et al., 2008 Fernandez et al., 2012), and 



the Complex C accretion model in Thorn et al. (2008). 

Based on Figures [I] and [2] and the discussion at the end of §2.2[ the ionized com- 
ponent of the HVCs extends beyond the immediate vicinity of the HI component. 
The calculation above includes only a factor of two on the HI mass for the warm 
ionized component directly around HVCs detected in Ha emission. Observations 
of ionized HVCs in absorption find an additional mass in more diffuse, warm ma- 
terial on the order of 10 8 M (Lehner & Howk, 2011). The ionized HVCs are at 
similar distances to the HI HVCs and also generally have a similar kinematic and 
spatial distribution, so to roughly account for the diffuse ionized gas the above 
best fit accretion rate can be scaled up by approximately a factor of two. Even 
with this additional ionized component, the accretion rate is still well below that 
needed by our Galaxy to maintain its current SFR over the next few Gyrs. Our 
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Galaxy may go through a low point in its gas content, and/or decrease its SFR 
until the arrival of the gas from the Magellanic System. 

For other galaxies, the amount of HI halo gas is generally between 10 7 ' 5-9 M , 
but this includes a combination of halo and disk-halo gas (within a few kpc). 
Some galaxies show very limited amounts of halo gas and substantial amounts 
at the disk-halo interface (e.g., Zschaechner et al., 2011 Heald et al., 2011). If 



accretion rates are estimated from the gas most akin to HVCs, Sancisi et al. (2008) ) 
calculates values of 0.1-0.2 M /yr for several galaxies from the HI component 
alone. Large covering fractions of ionized gas are detected through absorption 
line observations, and so more mass is certainly present in warm and warm-hot 
halo gas similar to the Milky Way. Including similar percentages of ionized gas 
as found for the Milky Way still leaves the accretion rates low compared to each 
galaxies SFR. The sources discussed below will provide additional fuel, but low 



redshift spiral galaxies may also gradually decrease their SFR's (e.g., Madau 
et al., 19981 iHopkins, 20041). 



6.2 Disk-Halo Interface 

Halo gas sources need to cool and integrate into the disk in order to become star 
formation fuel. Cooling and accretion at the disk-halo interface would be difficult 
to observe for galaxies beyond z ~ (i.e., no significant spatial and kinematic 
offset from the galaxy), and can be called 'quiet accretion'. This quiet accretion 
could explain observations at higher redshifts that show substantial SFRs and 



outflows with little, if any, evidence for accretion (e.g., Erb, 2008; Shapley, 2011 



Steidel et al., 2010). Cooling at the disk-halo interface also potentially explains 
the relatively constant amount of HI in the universe since z=3, despite the ongoing 
star formation 



(e.g., 



Prochaska & Wolfe, 2009 Putman et al., 2009a; Hopkins 



et al., 2008). Finally, quiet accretion may be linked to the fact that Hi-rich 



galaxies show quiescent disks (Wang et al., 2011). 



There is observational evidence for infalling structures at the disk-halo interface 
at low redshift. There are large HI features directly connected to galaxy disks in 
position-velocity space and the IVCs of the Milky Way have a clear bias towards 
negative velocities. In addition, the WIM layer of our Galaxy detected in Ha 



shows a net inflow as one examines the regions toward both poles (Haffner et al., 



2003 Putman et al., 2009b). A combination of infalling fuel and feedback can 



successfully model the lagging rotation with z-height found for numerous galaxies 



(Fraternali & Binney, 2008; Marinacci et al., 2011). The gradient of increasing 



temperature with increasing z-height (£ |2.2| &; £3.2) may also be a combination of 
the gradual cooling of halo gas, and the hot gas from the disk rising to larger 
z-heights than the cool gas. Given stellar metallicity constraints, the ratio of 
feedback material to fresh gas throughout time requires a full assessment. 

The interaction of inflowing clouds with existing halo gas is likely to be import- 



ant for cooling to occur at the disk-halo interface (Marinacci et al., 2010 Heitsch 



& Putman, 2009). As clouds move through the halo medium, they are gradually 



disrupted and slowed. In simulations of Heitsch & Putman (2009), the leftover 
warm over-densities of the HI clouds become buoyant as they approach the disk 
and are able to re-cool as they are compressed by the surrounding disk-halo me- 
dium. The small clouds that form in this process should be moving similarly to 
the gas surrounding it at this stage and may be the small cold clouds detected 
throughout the disk-halo interface of the Galaxy (q2~5J). When the clouds begin 
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to fall towards the disk again they will be moving through the surrounding diffuse 
medium and are likely to resemble the small warm intermediate velocity clouds 
that show net infall (Saul et al., 2012), or the filaments/ worms that emanate 
from the HI disk gas (e.g., Kang & Koo, 2007). 

The rate of accretion of different gaseous components at the disk-halo inter- 
face can be measured through a comparison of observations to high-resolution 
simulations. Preliminary examination of the MW simulation shows a significant 
amount of accretion of multiphase gas at this interface (Joung et al. (2012a); 
see also Murante et al. (2012[ )). This may indicate that quiet accretion at the 
disk-halo interface is an important fueling source for spiral galaxy disks. 



6.3 Disk Sources 



Obviously the gas in a galaxy's disk fuels its star formation, however there are 
two particular fuel reservoirs in the disk for which there is much left to learn. 
The first is the gas fed back into the ISM from evolved stars, and the second 
is outer disk gas that is transported to the inner star forming regions. AGB 
stars lose on the order of 30 — 50% of their mass back to the ISM according to 
stellar evolution models (Wachter et al., 2002). These evolved stars have been 
observed feeding star formation fuel back into the ISM with CO (e.g., Neri et al., 



1998[ |Castro-Carrizo et al., 2010) and HI observations (e.g., Gerard <fc Le Bertre 



2006 [ | Matthews fc Reid, 2007| | Putman et al., 2011a[ ). |Leitner fc Kravtsov (20lTf 



propose this type of stellar feedback may have provided a significant fraction of a 
spiral galaxy's fuel over the past few Gyrs. Galactic stellar metallicities and the 
deuterium abundance require the feedback material is supplemented with a source 
of lower metallicity gas. Understanding the methods with which the feedback 
material is fed back into the ISM, and the quantity, will help to determine the 
amount we need from external sources to maintain star formation. 

Another important question related to disk sources of fuel, is how material 
from gas-rich outer disks is transported to the inner star forming regions. A large 
reservoir of fuel is in the outer parts of the disk where few stars are currently 



forming ( Werk et al., 2010a; Thilker et al., 2007), and the flat metallicity gradients 



in gas disks indicate transport mechanisms are active (Werk et al., 2010b, 2011 



Bresolin et al., 2009 Spitoni fc Matteucci, 2011 ). The radial flow of gas in the 
disk would transport this fuel (Sellwood & Binney, 2002 Elson et al., 2011), as 



well as the potential return of material from a galaxy's warp to more central 
regions (Putman et al., 2009c Jozsa, 2011; Kawata et al., 2003). Perturbations 



by companions will result in mixing (e.g., Torrey et al., 2012) and can warp the 
gas disks (Binney, 1992). This will also create halo and disk-halo gas features 



and may be an important part of cycling accreted material inward. 



7 Future Directions 



We have begun to create maps of the multiphase gas in galaxy halos and de- 
termine its multiple origins and evolution. The relationship between halo gas 
and galaxy feedback and accretion mechanisms can be further clarified through 
additional studies in several areas. For instance, satellite gas is stripped into 
galaxy halos, but the percentage of halo gas that can be attributed to minor (or 
potentially past major) mergers remains to be determined. The IGM is rich with 
baryons, particularly near star forming galaxies, but our understanding of how 



27 



the IGM accretes onto galaxies is heavily based on simulation results. Feedback 
mechanisms impact galaxy halos, though how effective they are at transporting 
mass out of the disk and the primary mechanism responsible remains uncertain. 
Finally, although spiral galaxies have abundant fueling resources, the continuity 
of these resources and the physical mechanisms that lead to the conversion of the 
incoming fuel into molecular gas have yet to be fully determined. 

On the observational front, the next generation of telescopes and instruments 
will greatly aid our understanding of the flow of baryons in the Universe. Our 
knowledge of the link between the cold and warm gas detected in emission and 
absorption will become stronger with future HI maps with the Square Kilometer 
Array and the SKA pathfinders. Three-dimensional maps of the warm gas could 
be made with large area, high kinematic resolution integral field units (IFUs), 
and this would show the flow of warm gas relative to the disk. The hot halo gas 
remains an elusive component beyond the disk-halo region. Highly sensitive x-ray 
telescopes are needed for a direct detection, but in the meantime we can infer the 



presence of hot gas with broad warm-hot absorbers in the ultraviolet (e.g., Savage 



et al., 2011). The motion of the diffuse, hot halo component can potentially be 



calculated from the properties of the cold and warm clouds embedded within it. 
The characterization of the multi-phase halo gas in relation to a wide range of 
spiral galaxy properties (e.g., total and baryonic mass, SFR, environment) will 
provide significant insight into its origin. 

Stars in galaxy halos will be mapped with future large optical surveys, such 
as PanSTARRS and LSST. This will provide more information on the satellite 
accretion history and stars with a standard intrinsic color can be used for dust 
detection in the Galactic halo through reddening effects. Observing the halo stars 
for absorption lines from halo gas in the optical and ultraviolet will provide tighter 
distance constraints, and the lines can be used to examine the metallicity and 
ionization conditions of the halo gas. With additional estimates of the metallicity, 
halo gas origin scenarios can be further constrained (in particular the level of 
feedback), and the dust can be characterized. 

On the theoretical front, two of the most important steps are to continue to 
improve the resolution of the simulations and to study the effect of the many 
relevant physical processes, both separately and in combination, in a cosmolo- 
gical setting. Current cosmological models are unable to reach the length scale 
of many of the observed small HI clouds in the halo and at the disk-halo inter- 
face, nor can they accurately resolve the interaction of the multiple phases of 
gas with one another. These problems can be partially alleviated in the near 
future with models that are limited to a small region in space and time; however, 
then the relationship of the halo gas to galactic and cosmological processes on 
large scales is lost. Until supercomputers become fast enough to cover the large 
dynamic range required, clever combinations of cosmological and local models 
can be used to interpret the broad range of observed halo properties and dis- 
tinguish between various feedback and fueling mechanisms. Reproducing stellar 
metallicity distributions with different accretion scenarios will also be an import- 
ant feature of future simulations. The current chemical evolution models of our 



Galaxy are largely semi-analytic (although see Few et al., 2012; Kobayashi & 



Nakasato, 2011). Future state-of-the-art simulations can aim to track the metal 



flows in a galaxy, and through comparisons of observed and simulated stellar 
metallicities we can learn how feedback and accretion operate throughout time. 
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Fi gUXe 1: The distribution of HI (shaded clouds and plus symbols) and ionized high velocity gas 
(circles and diamonds) on the sky with color denoting the LSR velocity of the detection. The HI map 
was created with LAB data by |Westmeier (2007| l by removing the HI model of the Milky Way from 
Kalb erla fc Kerp (2009^. The plus s ymbols represent the the small HI compact HVCs detected with 
other data sets (|Putman et al., 2002| |de Heij et al., 2002| >, open circles are the O VI absorption line 
detections jSembach et al,, 2003| l, solid circles are the Si absorption line detectio ns ( |Shull et al., 2009"t , 
and diamonds are the ionized HVCs (IHVCs) at J;15 kpc I Lehner & Howk, 2011 1. Multiple colors for a 
symbol indicate multiple absorbers along the line of sight. The positions of the major HI HVC complexes 
(with S to indicate the Smith cloud) and the LMC, SMC, M31 and M33 galaxies are noted. The AC 
Shell appears in this map, but is considered an intermediate velocity cloud. 
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Figure 2: The distribution of HI (shaded clouds and plus symbols) and ionized high velocity gas 
(circles and diamonds) on the sky with color denoting the GSR velocity of the detection. Otherwise this 
figure is the same as Figure [I] with the exception of showing the view from the North (left) and South 
(right) Galactic Poles at the bottom. 
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Figure 3: A 3D, interactive model of the Galactic HVC complexes. The Galactic disk is shown with a 
25 kpc radius and a 1 kpc thickness. The central regions of our Galaxy are indicated with a gray 3 kpc 
radius sphere, and the Solar position is shown with a yellow sphere. HVC complexes are placed at the 
average of their most stringent distance constraints (j ]2.1.1| . The Magellanic Stream is shown starting 
at a distance of 50 kpc, and increasing linearly to 150 kpc at a Magellanic longitude of 150 degrees, 
consistent with the models of |Besla et al. (20 10). The positions of the LMC and SMC are noted with 
gray spheres. Cloud thickness is an arbitrary function of column density along the line of sight. If the 
Figure does not appear interactive, please view with a current version of Adobe Reader and click and 
drag to rotate. 
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Figure 4: Left: Head-tail HI halo clouds with the outer contour at N H i = 2 X 10 18 cm~ 2 (15.5' res- 
olution; |Putma^^t^L7^011b]| . Right: A HI HVC that shows signatures of interaction with the halo 
medium from the edge features (4' resolution; |Peek et al., 2007^ . 
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Figure 5: A cartoon representing the disk-halo interface given the observations discussed in [2.5 



and 



*j3] The temperature gradient for the diffuse gas is shown with the pink color bar, while the kinematic 
gradient of the denser neutral gas (HI) is shown by the blue-gray shading. There are distinct populations 
of large IVCs and small discrete clouds, but in reality the edges of the HI features are not this sharp. 
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Figure 6: Comparisons of HI halo gas for M31 (left) and NGC 891 (right) with that in the Milky Way 
if viewed with the same inclination. The exterior MW HI view is reconstructed using the distances to the 
13 known HVC complexes (white) and assuming a distance for the MS/LA complex (red) ranging from 
50 to 150 kpc using the models of ( |Besla et al., 2010| l. The MW disk (blue), including the disk column 
density, flare, and warp, is reconstructed from the work of |Levine et al. (2006}. Contours are set to m atch 
the observations: 0.1, 0.3, 1, 3.2, 10, 32, 100, 316, 1000 x IP 1 " cm " 2 for the|Westmeier et al. (2007} M31 
data and 0.5, 1, 2, 5, 10, 20, 50, 100, 200, 500 x 10 19 cm" 2 for the |Oosterloo et al. (2007} NGC 891 data, 
respectively. The large cloud above M31's disk is unlikely to be associated with M31 ( |Westmeier et al.,| 
2007}, and see |Ibata et al. (2007] and |Mouhcine et al. (2010} for stellar features around M31 and NGC 
891, respectively. Viewed externally, the Magellanic system is clearly the dominant feature in our HI 
halo and our HI disk is comparable to M31 and NGC 891. We find that, without velocity information, 
the non-Magellanic HVCs would be hard to distinguish from the flared, warped disk. 
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Figure 7: A cosmological AMR simulation of a Milky Way galaxy at z=0.3 (Joung et al., 2012a I. 
We chose t = 10.3 Gyr as a relatively low redshift time slice, as it clearly shows examples of three 
different origins of neutral gas in the halo as indicated by the numbered boxes in the top panel. Box 
1 represents the interaction between outflowing gas and accreting halo gas, box 2 shows a filamentary 
cool flow feeding the disk, and box 3 represents the accretion of three satellites. The cut-outs show the 
total gas density (N), temperature (T), stellar density (£*) and metallicity (Z). The features become 
less distinct in the total gas density plots due to the large amount of multiphase gas in the halo. 
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Fi gUrG 8: Top: A diagram based on the MW simulation at z=Q to a radius of 100 kpc where the 
top hemisphere shows the outflow due to SN-driven winds (pink) and the bottom shows the inflow of 
cooling (green to blue) and heating (green to red) filamentary flows and a satellite (black). Bottom: The 
density-temperature diagram for cells in the MW simulation at z=0 with color coding corresponding to 
the top panel and main features noted. The directions of increasing pressure and entropy are shown with 
the arrows in the upper right. The cells were randomly selected from the simulation with the probability 
of selection proportional to gas density. 
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Figure 9: The accretion of HVC complex mass with time for the fit (dark shades) and maximal (light 
shades) accretion velocities vq. Dashed lines indicate accretion rates of 0.08 M^/yr and 0.4 Mq/yt, 
respectively. Several specific HVC complexes are noted for reference. The gap at recent times may be 
due to selection effects for HVCs close to the disk. 



